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ABSTRACT 



Fluctuations in the horizontal component of the earth's 
magnetic field were measured at a land site (La Mesa) and 
on the floor of Monterey Bay at a depth of 70 meters. The 
measurements span a 19-hour period (August 17/18, 1982). 
Assuming a circular polarization basis set, analysis was 
conducted to determine the coherence between the two sites 
in the frequency range .05-10 Hz. Preliminary analysis 
indicates a high degree of coherence from .05-. 4 Hz. Also, 
a study of the polarization of the micropulsation field at 
each site was conducted in terms of Stokes parameters. 
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I. INTRODUCTION 



This thesis research is part of an ongoing effort being 
conducted by the Naval Postgraduate School to obtain a long- 
term data base of ELF electromagnetic noise on the ocean 
floor. The result of this project should be a better under- 
standing and interpretation of that noise and its sources. 

Magnetic noise is important to the geophysicist and to 
the Navy. Magnetic techniques are applied in exploration 
for minerals , the prediction of communications propagation 
and to the modeling of the sphere on which we live. Geo- 
magnetic studies are of interest to the Navy for their 
application to submarine detection, submarine communications 
and to mine warfare. 

The objective of this thesis within the broader context 
of the overall project was to design and validate a computer 
code that would facilitate data manipulation and analysis. 

To that end, the present computer code calculates for two 
sites (°ne on the ocean floor, the other on land), the 
individual site's power spectral densities and Stokes 
parameters 0,1, 2 and 3 for two of the individual site's 
orthogonal components in the horizontal plane. The code 
also computes the coherence of these two orthogonal compo- 
nents and their phase difference. Lastly, the program 
computes the right and left circular polarization power 
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spectral densities of the sites and the coherence between 
the seafloor and land-based sites. 

As discussed, the present code applies only to two 
sites and two orthogonal geomagnetic components at each 
site. However, extension to several sites and/or components 
should be an easy undertaking due to the style of the present 
program. 
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II. BACKGROUND 



A. THE MAIN GEOMAGNETIC FIELD 

Historical observation and modern measurement demonstrate 
that the principal source of the geomagnetic field (approxi- 
mately 50,000 nanoteslas) is beneath the surface of the 
earth. An easy- to-picture causal interpretation of this 
phenomenon is painted by imagining the earth with a dipole 
magnet near its center as pictured in Figure 2.1. The 
magnet axis is tilted at an angle of 11.5 degrees from the 
earth's rotation axis. This model provides a prediction 
accuracy to a few percent of the measured geomagnetic field. 
Figures accurate to less than one percent of the earth's 
measured field can be obtained using more sophisticated 
models [Ref . 1 ] . 

The generally accepted physical process for the genera- 
tion of the geomagnetic field is that of the self-sustaining 
dynamo [Ref. 2 ]. Figure 2.2 presents a simple picture of 
such a dynamo. The physical requirements of this theory 
are mostly met in the interior of the earth. For example, 
the earth's center is a large conducting system. Addition- 
ally, we can assume that the earth's molten core is of a 
differential rotational nature with turbulent upwelling. 

Energy sources for such a dynamo are most likely avail- 
able in the earth's core due to a number of processes. This 
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Figure 2.1 Eccentric Dipole Model of Geomagnetic Field 
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Figure 2.2 Simple Disk Dynamo 

all boils down to mean that the motion of a conductor, such 
as the earth's molten core, in a magnetic field produces a 
current which in turn produces a magnetic field in support 
of the original magnetic field [Ref. 3] . 

A simplified picture of the dynamo theory originates 
by assuming that in the distant past a magnetic "seed" field 
was captured by the earth's conductive core. The differen- 
tial rotation of this core "winds up" the field into an 
intense azimuthal field, such as would happen to a twisted 
rubber band. Turbulent upwelling of the core tends to 
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carry this azimuthal field outward, converting it to a 
toroidal helical field within an associated azimuthal 
circulating current. This current in turn generates the 
earth's dipolar field. 

An additional advantage of the self-sustaining dynamo 
theory lies in its ability to explain the well -documented 
reversals of the geomagnetic field. It has been postulated 
that two or more coupled "dynamos" as shown in Figure 2.3 
can cause reversals of the main geomagnetic field over 
geologic time scales [Ref. 4] . 




Figure 2.3 Twin Disk Dynamo 
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B. GEOMAGNETIC FLUCTUATIONS 



The geomagnetic field changes with time. The changes in 
the field have periods ranging from fractions of seconds to 
thousands of years. The very slow variations (periods of 
years or more) are not of interest here but are known as 
secular variations and are thought to be geologic in origin. 

The geomagnetic field changes of interest here are those 
fluctuations that occur over the period of a day or less. 

It is customary to view the magnetic field as the stable 
main field with these comparatively rapid fluctuations 
superimposed on it. Figure 2.4 presents the spectrum of 
the superimposed fluctuations. 

The fluctuations (Figure 2.4) that are superimposed on 
the main field range in intensity from .05 to 50 nanoteslas 
and have various causes. Some major contributors to these 
fluctuations are in the ionospheric currents overhead and 
in magnetic storms which probably originate in solar 
eruptions . 

Finally, arriving at the point of interest here, are 
those disturbances to the earth's magnetic field that 
propagate through the earth's magnetosphere as hydrodynamic 
waves. These small perturbations cover a wide range in 
amplitude and duration and are known as micropulsations. 

Figure 2.5 presents a graphical display of these disturbances. 
Micropulsations are thought to be the result of fluctuations 
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POWER DENSITY (nT) 




Frequency (Hz) 

Figure 2.4 Power Spectrum of Geomagnetic Disturbances 
Observed on the Surface of the Earth. 



[Claudis , Davidson, and Newkirk, 1971] 
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100 .01 

Period (sec) 

Figure 2.5 Field Strength of Micropulsations 
[Campbell, 1966] 

in the incident solar wind affecting wave -particle interaction 
in the magnetosphere [Ref. 5] . As they propagate through 
the conductive layers of the atmosphere, they are trans- 
formed into electromagnetic waves via ionospheric currents, 
which we detect at the earth's surface as magnetic field 
fluctuations. There are two major categories of micro- 
pulsations: continuous and irregular. 

Smoothly varying, periodic pulsations are classified as 
(Pc) micropulsations. Their distinct periodicity is indica- 
tive of a resonance inside the magnetospheric cavity. One 
model proposed to explain this occurrence of smooth, periodic 
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pulses is that the resonance is a result of an oscillatory 
hydromagnetic wave. This wave is associated with groups 
of ions which "bounce" back and forth in the magnetosphere 
[Ref. 6] . 

These continuous micropulsations are further classified 
according to their periods as follows: 

1. Pci: [0*15-5 second period) 

Known as "Pearls," these micropulsations 
are generated by the cyclotron instability 
of the energetic protons. They have been posi- 
tively correlated with solar disturbances and 
occur during daylight hours in the auroral 
zone and during night and early morning hours in 
the midlatitudes. Their average amplitude is 
approximately 0.05-0.1 nanotesla [Ref. 7]. 

2. Pc2: (5-10 second period) 

This is a diurnal phenomenon that shows some 
positive correlation with solar activity and the 
seasons. They usually decrease in their period 
as magnetic activity increases. Their average 
amplitude is 0.1-1 nanotesla. 

3. Pc3: (10-45 second period) 

This micropulsation is similar to Pc2 pulsations 
except for the difference in period. 

4. Pc4: (45-150 second period) 

Sunspot activity appears to have an effect on 
Pc4 pulsations. The Pc4 pulsations also change 
their period seasonally with an average amplitude 
of 5-10 nanotesla. 

5. Pc5 : (150-600 second period) 

These large-scale pulsations typically occur 
during morning and evening with amplitudes of 
10-100 nanotesla. Their duration shows a strong 
geomagnetic latitude dependence. 

Irregular (Pi) micropulsations appear to be associated 

with polar magnetic disturbances and are observed during 

night-time hours. The theory is that these micropulsations 
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arise from disturbances of the ionosphere due to particle 
bombardment from the magnetosphere. Increased solar radia- 
tion levels also appear to have an effect on their existence. 
Irregular micropulsations are classified as follows: 

1. Pil : (1-40 second period) 

These pulsations have a small amplitude and 
usually occur at night and early morning. 

They vary in intensity from 0.01-0.1 nanotesla. 

They demonstrate a positive correlation with 
auroral disturbances. 

2. Pi2 : (40-150 second period) 

Pi2 pulsation amplitude range is on the order 
of 1-5 nanotesla. They usually occur during 
the early evening hours but may continue through- 
out the night. The periods of these pulsations 
decreases with increasing magnetic activity. 

The Pc5 and Pc4 micropulsations last several hours. 
However, the Pc3-1 micropulsations and irregular pulsations 
have a maximum duration of approximately one hour but may 
cease after only a few minutes [Ref. 8]. Above 3 Hz, the 
most pronounced field variations are the earth ionosphere 
cavity resonances of Schumann resonances [Ref. 9] . 



C. OCEAN WAVE EFFECTS 

In addition to all of the geomagnetic fluctuations 
discussed above, if working in the ocean environment, one 
must take into account hydrodynamically induced magnetic 
fluctuations. These fluctuations are internally produced 
resulting from the movement of a conductor (seawater) in 
an external magnetic field. The effects of waves have been 
calculated [Ref. 10] and are pictured in Figure 2.6. 
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DEPTH IN METERS 




Figure 2.6 Induced Magnetic Field per Meter Amplitude 
of the Surface Wave 

[Weaver, 1965] 
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Weaver's equations have been modified to incorporate bottom 
reflection and transmission with the result that the effect 
of bottom decreased the magnitude of the induced magnetic 
field by less than 3 dB for frequencies greater than 0.16 Hz 
in water depths less than 500 meters [Ref. 11]. 
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III. COLLECTION SYSTEM 



A. DATA COLLECTION SYSTEM 

The data acquisition system is illustrated in Figure 3.1. 
Major components consist of: 

Cl) coil antenna sensors (2) 

(2) preamplifiers (2) 

(3) signal conditioners (2) 

(4) pulse code modulation system (1) 

(5) radio transmitter (1) 

(6) radio receiver (1) 

(7) tape recorder (1) 

For the work reported in this paper, two systems consisting 
of components (1) through (5) were deployed. They were 
similar in most respects with the major exception that the 
ocean equipment was enclosed in watertight glass spheres. 

A thorough description of the data acquisition system is 
given in Reference 12. 

1 . Coil Antenna Sensors 

Each coil antenna sensor consists of two continu- 
ously wound coils which are constructed of about 5460 turns 
of 18-gauge copper magnet wire. Glass spheres are used to 
maintain the watertight integrity of the coils and all other 
immersed equipment for the ocean configuration. See Figure 
3.2 for the ocean system. The coils are mounted orthogonally 
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Figure 3.1 Data Collection System 








28 




in a nonmagnetic frame with each coil having a resistance 
of 120 ohms and an inductance of approximately 9.31 Henries. 

2 . Preamplifiers 

The preamplifiers are model 13-10A low-noise ELF 
amplifiers manufactured by Dr. Alan Philips. The final 
stage of the amplifier contains an active low-pass filter 
which provides a sharp cutoff for frequencies above 20 Hz. 
Frequencies above 10 Hz are not examined in this paper. The 
overall amplifier gain for inputs of less than 2.5 milli- 
volts is 60 dB in power. The purpose of the filter is to 
remove 60 Hz noise and to prefilter the data before 
digitization . 

3 . Signal Conditioners 

The signal conditioners receive analog signals from 
the preamplifiers, amplify them by about 30 dB and limit 
signals with peak amplitudes of 7.5 volts (modifiable) from 
entering the pulse code modulation system. 

4 . Pulse Code Modulation System 

The pulse code modulation system was designed and 
manufactured by Dr. Robert Lowe, Lowecom, Inc. The system 
offers a 15-channel analog input capability with selectable 
sampling rates of 2, 4, 8, 16, 32, 64 and 128 samples per 
second. The PCM system uses a crystal oscillator and inte- 
grated circuitry to develop the clocking pulses, a 16-channel 
analog multiplexer, a 16-channel, 12-bit analog-to-digital 
converter as the basis for the pulse coding. The crystal 
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clock oscillator loses 1 bit in 1 million. The clock pulse 
gates the analog multiplexer, the analog-to-digital converter 
and additional circuitry that form the pulse code words. The 
PCM output is organized in frames with each frame headed by 
a synchronization word which is followed by the pulse-coded 
samples from PCM channels 1 through 15. 

5 . Transmission and Recording 

After encoding, the data is transmitted via VHF 
radio to a Naval Postgraduate School located receiver and 
recorded on an instrumentation tape recorder. For the data 
of this paper, a sample rate of 32 samples per second was used 

B. DATA ANALYSIS EQUIPMENT 

1 . Analysis Hardware 

The PCM data is fed into a PCM decoder which places 
the decoded PCM onto a 9- track 800 BPI computer tape. This 
computer tape is read into the Naval Postgraduate School IBM 
3033 computer for manipulation and analysis. 

2 . Data Analysis Software 

The IBM 3033 is used to convert the data to usable 
form and then to perform any operations on the data as desired 
At present the computer code is written in Fortran IV and 
converts the digital data to voltage, Fast Fourier trans- 
forms the data, calculates power spectral densities, deter- 
mines Stokes parameters 0 through 3 and provides phase and 
coherence figures. Appendix A contains test verification 
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of the code. Appendix B contains the code that prepares 
data from one site for storage in the IBM 3033 Mass Storage 
System for later recall and comparison with another site. 
Appendix C contains the code that prepares and manipulates 
the code from a second site, retrieves from the Mass Storage 
System that data previously prepared and manipulates the 
two data sets simultaneously. 
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IV. COMPUTER CODE 



A. INTRODUCTION 

As mentioned earlier, one of the intentions of this 
effort was to develop a code that would permit the comparison 
and/or simultaneous geomagnetic analysis of two different 
sites. The present data collection system provides for the 
recording of transmitted data from those sites onto different 
tracks of the instrumentation tape recorder. Since only one 
track of that analog tape can be placed on one digital 
computer tape, a means of storing one channel had to be 
devised. As can be seen by a quick review of Appendices B 
and C, the present technique that allows this simultaneous 
manipulation of data from two sites is to read the computer 
tape from one site and to process that data to the point 
where it is ready for analysis. These results are then stored 
in the IBM 3033. Subsequently, the same routine (with modifi- 
cations to fit the parameters of the second site) are applied 
to data read from the second site computer tape. Now both 
sets of data are in the same form and ready for manipulation. 

B. DATA ANALYSIS SOFTWARE 

The two programs (one for each of two sites) are 
essentially the same with the exception that one program, 
program B, is a short version of the other, program A. Thus, 
the next paragraphs will discuss in detail the longer program. 
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This description fits equally well to program B except that 
program B ends after applying the system transfer functions 
and placing these results in the Mass Storage System. 

1 . The Main Program 

The main program is divided into sections which 
perform the following functions: 

* Data input 

* Fourier analysis of time series data 

* Application of system transfer function. Note: At 

this point, program B reads the data into the Mass 

Storage System and ends. Program A, at this point, 

retrieves that same data from Mass Storage. 

* Data averaging 

* Calculation and plotting of: 

(1) power spectral densities of the individual 
site geomagnetic orthogonal components. 

(2) orthogonal component coherence of the geo- 
magnetic field of one site. 

(3) Stokes parameters 0 through 3 of the individual 
site orthogonal components. 

(4) single site right and left circular polarization 
power spectral densities. 

(5) single plane coherence of right and left 
circular polarization between sites. 

a. Data Input 

The data is read into the IBM 3033 via the sub- 
routine '‘read" which was written by Tim Stanton of the 
Naval Postgraduate School. The "read" subroutine reads the 
PCM data from the computer tape and converts it into an 
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integer from 0 to 4095. The main program than converts 
this integer to voltage. See Appendix B or C. 

b. Fourier Analysis 

The next section of the code utilizes a sub- 
routine developed at the Naval Postgraduate School and 
called "fourt." "Fourt" Fast Fourier transforms the voltage 
data from the time domain to the frequency domain. See 
Appendix A for validation of this routine . 

c. Transfer Functions 

After the Fast Fourier transform, the program 
applies the transfer functions for the system. Development 
of the transfer functions is contained in Reference 12. 

The transfer functions are applied up to 10 Hz as straight 
line functions. The data is now in nanoteslas. 

d. Power Spectral Density 

This quantity is computed by multiplying the 
previously calculated spectral density summation by the 
observation period and dividing by the number of data 
groups (typically 20) that the computer tape was segmented 
into, providing an average power. 

e. Coherence 

The next steps of the program determined for one 
site the geomagnetic orthogonal components Stokes parameters 
0, 1, 2, and 3, and the phase of the components and their 
coherence. Next, right and left circular polarization 
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parameters are computed for the individual sites followed 
by calculation of the coherence of these parameters between 
the two sites. 



C . PARAMETER THEORY 

1 . Stokes Parameters 

A useful tool for determining the state of polariza- 
tion of a wave are the Stokes parameters 0, 1, 2, and 3 
[Ref. 13]. 

Two linearly polarized plane waves with field 
vectors in the directions of and can be described as: 



. -V 



£ = p F P (i k>x " iwt ) 

*1 £ 1 h l e 



E 2 - e 2 E 2 e 



(i %‘X - iwt) 



These waves can be combined to give a general plane wave 
propagating in the direction $ = kn: 



£ (x,t) 




+ 



/s 






1 *" 

k *x 



iwt) 



CD 



If E 1 and E^ have the same phase, (1) represents a linearly 
polarized wave with its polarization vector making an angle 

0 = tan" (E /E^) with E^ and with a magnitude E = *4;^ + E^ . 
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If E^ and have different phases, the wave (1) is 
elliptically polarized. 

In the simplest case, circular polarization, E^ and 
E 2 have the same amplitude but differ in phase by 90°. The 
wave (1) is then: 

I (5,t) - E 0 ± e 2 J e ^ • iut > (2) 

Then the components of (2) are: 

E (x,t) = E cos (kz - u t) 

X o 

( 3) 

E (X,t) = E q sin (kz - tot) 

Thus, at a fixed point in space, the fields (3) are swept 
around circularly at the frequency to. For (e^ + , the 

rotation is counterclockwise if the observer is facing the 
oncoming wave. This wave is described as having positive 
helicity. (e^ - ie 2 ) has negative helicity or is rotating 
clockwise when looking into the wave. 

The two circularly polarized waves (2) form an 
equally acceptable set of basis fields for description of 
a general state of polarization by introducing the complex 
orthogonal unit vectors : 
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Then a general representation of (1) is: 






S (J,t) = CE + * E. h e (i k ' x 



iu)t) 



(4) 



where E + and E_ are complex amplitudes. If E + and E have 
different amplitudes but the same phase, (4) is an elliptically 
polarized wave with principal axes of the ellipse in the plane 
defined by and e ^ . 

Thus, the polarization content of a plane wave is 
known if it can be written in the form (1) or (4) with known 
coefficients (E E 2 ) or (E + , E_) . Conversely, how can we 
determine from observations on a wave the state of its 
polarization? The answer lies in the use of the four Stokes 
parameters. Their measurement or calculation determines 
the state of polarization of a wave. 

For a wave propagating in the z direction, the scalar 
products : 
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are the amplitudes of radiation respectively, with linear 
polarization in the x direction, linear polarization in the 
y direction, positive and negative helicity. The squares of 
these amplitudes provide a measure of the intensity of each 
type of polarization. Phase information is obtained from 
cross products. If we assume the scalar coefficients of 
(1) and (4) to be a magnitude times a phase factor: 



and use the circular polarization basis (e + , e_) , we have: 




E 




E 



a. e 



E = a e 




=2Re[(e + • E + ) (e_ • E) ] - 2a + a_ cos(<5_ - <S + ) 



S 2 = 2Im[(e + ♦ E + ) * (e_ • EJ] = 2a + a_ sin(6_ - 5 + ) 
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2. Coherence 



If a plane wave of frequency v is propagating in the 
z direction, it has components: 

B x (t) - a x (t) e i< ^ 1 ' 12lrvt 

Ey (t) = - i2,tvt 

that are mutually orthogonal and at right angles to the 
direction of propagation [Ref. 13]. Then the complex 
degree of coherence can be written as : 



U 



xy 



J 



/T 



xx 



xx_ 

■ /J 



yy 



where J ,... are the elements of the coherency matrix 

XX 
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E E 
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E E 






i+i-if* 


X X 


x y 
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E E 


* 

E E 
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i 

i 

X 


y y 




1 2 



The complex factor u is the complex degree of coherence and 

xy 

is a measure of the correlation between components of the 
field vector in the x and y directions. 
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D. SOFTWARE TESTING 



The basis set for the code calculations is that of a 
circularly polarized wave. Thus a circularly polarized 
wave was chosen as the test input. The wave was created by 
generating a 4 Hz sine wave for sites A and B x coil input 
and 4 Hz cosine wave for sites A and B y coil input. The 
results of this test are valid only at the input frequency 
as a measure of the code veracity. One of the characteris- 
tics of any FFT routine is that small residuals are created 
throughout the spectrum of interest. These residuals are 
primarily caused by the finite observation time. These 
small-valued artifacts will be carried throughout all subse- 
quent parameter calculations and will appear in the test 
plots. This inclusion of unexpected values at frequencies 
other than the test input of 4 Hz should be considered as 
the result of the artificiality of the single value test 
input. Input of real world data that covers a full spectrum 
of values obscures the "FOURT" routine propensity to create 
these artifacts. 

1. Sites A and B Orthogonal Component 
Power Spectral Density Test 

Spectral density is defined as: 



S 



x 




(x^Cx.)* 



for site A x coil where x is the Fourier Transform of some 
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input to the x coil. It follows that the average power 
spectral density is: 



PSD o 
xa 



Time 

~N 



n 

S 

1 






S xa ' Time 



where time is the observation time. Thus, for the test 
inputs of sine 4 Hz into the x coils at sites A and B and 
the cosine 4 Hz into the y coils, we expect and obtain (see 
Figures A.l, A. 2, A. 11, A. 12) a peak magnitude of .5 at 
4 Hz at each site. 

2 . Sites A and B Orthogonal Component Coherence Test 
Coherence as defined above is : 



U 



xya 



S 

xya 



/3 



xa 



/5 — 
ya 



for two 



orthogonal components where 



S 

xya 



is 



defined by: 



S 

xya 



Time 

N 



n 

Z 

1 






With X and Y our transformed inputs and considering the same 
inputs, we expect to see a coherence of magnitude 1 at a 
frequency of 4 Hz and its harmonics. Figures A. 8 and A. 18 
demonstrate the validity of the coherence calculations of 
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the code. Also the value of 1 for the test case coherence 
at frequencies other than 4 Hz is a manifestation of the 
residual effect discussed above. 

3 . Sites A and B Stokes Parameter 1 Test 

Using the definitions developed in preceding 
paragraphs, Stokes 1 is: 



S la = 2a + a . sin ( 6 _ ‘ <$+) " 2[S x - S v ] = 0 



x y- 



for S = S as is the case for this test input and as is 
x y 

shown in Figures A. 4 and A. 13. 

4 . Sites A and B Stokes Parameter 2 Test 

Again using definitions developed in the preceding 
paragraphs, Stokes parameter 2 is: 



S ?a = 2a + a_sin (6_ - 6 + ) = 0 for (6_ - 6 + ) = 0 



Figures A. 5 and A. 15 validate the code for this calculation 

5 . Sites A and B Stokes Parameter 3 Test 
Stokes 3 is : 



S 3a = a + - a ! = - ( • 707) 2 - ( . 707) 2 = -1 



Figures A. 6 and A. 16 validate the code for this calculation 
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6. Sites A and B Stokes Parameter 0 Test 



Stokes 0 is : 



S_ = a? + a 2 = 2 [S „ = S . ] = 2 
oa + - 1 xa xb J 



(See Figures A. 3 and A. 13.) 

7 . Sites A and B Phase Difference of x and y Test 

The phase difference of sine 4 Hz versus cosine 4 Hz 
is 90° or 1.57 radians. (See Figures A. 7 and A. 17.) 

8. Sites A and B Two Orthogonal Component Power 
Spectral Density with Right and Left Circular 
Polarization Basis - Set Test 

The development of the spectral densities for this 
case is similar to the development of spectral density for 
a single component. Starting with the cross spectra density 
definition for right and left circular polarization 



BP 



xya 



(X + iY) (right polarized) 



BM 



xya 



(X - iY) (left polarized) 



It follows that the spectral density is: 




(right polarized) 




(left polarized) 
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The power spectral densities will be: 



SP 



xya 



Z SP 

N ^ xya 



SM 



xya 



liB®. ? SM 
N ^ xya 



where the time is the observation time. Examples of this 

calculation with a test input of 4 Hz right circularly 

polarized wave are in Figures A. 9, A. 10, A. 19 and A. 20. 

9. Coherence of Sites A and B Right and Left Circular 
Polarization Basis Set Test 

Right and left circular polarizations are: 



SP 



COP = 



abxy 



/3P7 /SET 
a b 



(right polarization) 



SM 



COM = 



abxy 



/SM~ /5MT 

a b 



(left polarization) 



where SP , and SM , are the cross spectra densities of 
abxy abxy 

the two orthogonal components at the sites A and B. 



SP = IAE!£. £ (BP ) (BP ) (right polarization) 

abxy N 1 xya' xya 7 
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(left polarization) 



SM abxy * T 1 \ ( BM xya> < BM xyb> 

The test case of sine 4 Hz as the xa and xb inputs versus 
cosine 4 Hz as the ya and yb inputs yields a coherence of 0 
for left circular polarization and 1 for right circular 
polarization. (See Figures A. 20 and A. 21.) 
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V. ANALYSIS 



Preliminary, "eyeball" analysis was conducted on data 
collected over a period from 1122 local 17 Aug 82 to 0639 
local 18 Aug 82. The timeframes inspected within this 
window were : 



TIME (local) 




DATE 




1122-1254 


17 


Aug 


82 


1301-1408 


17 


Aug 


82 


2028-2200 


17 


Aug 


82 


0121-0251 


18 


Aug 


82 


0259-0431 


18 


Aug 


82 


0507-0639 


18 


Aug 


82 



Original intentions were to collect data over a full 24-hour 
period. However, the batteries at the Monterey Bay site went 
dead after 20 hours. Of the 20 hours of collected data, 
eight hours appeared dominated by system noise and were not 
included in the analysis. However, the remaining data of 
approximately eight hours spans a 24-hour period and provides 
evidence of some obvious trends that lend themselves to the 
cursory analysis conducted here. Representative plots of 
the calculated parameters are in Appendix D. Technical 
Report NPS 61-83-005 contains a complete set of plots of the 
calculations that span the full 20 hours of data collection. 

This complete set of plots was used in the analysis that follows. 
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A. POWER SPECTRAL DENSITIES 



Both sites exhibit a marked decrease in power spectral 
density with frequency. There is an obvious exception to 
this observation at 4 Hz at the land site. The sharp peak 
at 4 Hz is attributed to 60 Hz aliasing which is a result 
of the Fast Fourier Transform (FFT) subroutine. Fast 
Fourier algorithm construction is such that data being 
transformed may be convolved back onto itself [Ref. 15]. 

This convolution is routinely avoided by using a sample 
rate that is twice the upper frequency of the domain of 
interest. A sample rate of 32 samples per second was used 
in this experiment and should have precluded aliasing of 
the frequency of 60 Hz, especially when combined with the 
amplifier filter cutoff of 20 Hz. However, the 60 Hz 
component of the magnetic field was measured and found to 
be 4 nanotesla/Hz which is 10 s higher than magnetic noise 
data [Ref. 16]. The 20 Hz filter cannot completely filter 
out signals of this magnitude. Thus the 4 Hz signal is a 
result of aliasing of the strong ambient 60 Hz signal. 

B. STOKES PARAMETERS 

All Stokes parameters demonstrate reasonable agreement. 
That is, Stokes 0 is obviously a 2 + a* and Stokes 3 is 
obviously a 2 - a 2 as shown in Figures D.3, D.6, D.13 and D.16. 
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1. Ocean Site 



The ocean site Stokes 2 and 3 parameters exhibit a 
time dependence. That is, in the early a.m. hours the large 
magnitude of Stokes 2 is the noticeable feature. This 
magnitude decreases at sunrise and through the day until it 
reasserts itself after sunset. A review of Stokes parameters 
2 and 3 over the total observation time (see Technical Report 
NPS-61-83-005) indicates that this falloff is independent of 
frequency. The increase in magnitude of Stokes 2 after 
sunset is also independent of frequency with the exception 
of some very discrete spectral components. 

2 . Land Site 

The land site Stokes 2 and 3 parameters show a time 
dependence with a shorter cycle. That is, the magnitude of 
Stokes 2 is large near sunset, followed by its decrease 
during the morning hours. The early afternoon hours reflect 
an increase in Stokes 2 until after sunset. Figures 5.1 
through 5.4 show the decreased magnitude of Stokes parameters 
2 and 3 for the morning hours. 

C. X AND Y COIL PHASE DIFFERENCE 
1. Ocean Site 

The use of Stokes parameters in the calculation of 
the phase difference dictates that changes in the time 
series averaged Stokes parameters should be evident in the 
x and y component phase difference. This is the case as 
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Figure 5.1 Stokes Parameter 2, La Mesa, 1122-1254 Local, 17 Aug. 82 
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Figure 5.2 Stokes Parameter 3, La Mesa, 1122-1254 Local, 17 Aug. 82 
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Figure 5.3 Stokes Parameter 2, La Mesa, 0507-0639 Local, 18 Aug. 82 
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Figure 5.4 Stokes Parameter 3, La Mesa, 0507-0639 Local, 18 Aug. 82 



shown in Figures D.7 and D.17. The phase difference from 
sunset to sunrise shows a definite decrease with an "eye- 
ball" average of .2 to .4 radians. Daylight hours show 
a difference ranging from .5 to .3 radians with the greatest 
difference near local noon. 

2 . Land Site 

The land site phase difference reflects its Stokes 
parameters with a difference of about .3 radians in the 
higher frequencies at all hours and a difference of .1 
radians during the daylight hours. 

D. X AND Y COIL COHERENCE 

1 . Ocean Site 

Time dependence is again a strong factor with high 
degrees of coherence in the ocean site x and y coils being 
greatest during hours of darkness and diminishing at sunrise. 
Figures 5.5 and 5.6 are examples of this phenomenon. The 
exception to this observation is at local noon when the 
coherence of the ocean coils is again very high but diminishes 
before and after the noon time data. 

2 . Land Site 

The land site x and y coil coherence reflects the 
previous parameter cycle time dependence. There is a high 
degree of coherence in the lower frequencies at sunrise 
that is followed by diminishing coherence in the higher 
frequencies. Local noon shows an improved coherence at 
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figure 5.5 Coherence of x and y Coils, Mtry Bay, 18 Aug. 82, 0121-0251 Local 
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Figure 5.6 Coherence of x and y Coils, Mtry Bay, 18 Aug. 82, 1301-1408 Local 



these frequencies but is followed by a degraded coherence 
after sunset. 

E. CIRCULAR POLARIZATION POWER SPECTRAL DENSITY 
There is no marked contrast between right and left 

circular polarization power spectral densities. Both 
calculations provide similar results. 

F. COHERENCE OF CIRCULAR POLARIZATION 

No preference for right versus left circular polarization 
is discernable . The most obvious trend is that of a high 
degree of coherence in the frequency range .05 Hz to .4 Hz 
as shown in Figures 5.7 and 5.8. Generally, coherence of 
the two sites decreases thereafter with frequency. 
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FREQUENCY (IIZ) 

Figure 5.7 Coherence of Right Circular Polarization Monterey Bay/La Mesa 
1122-1254 Local, 17 Aug. 82 
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Figure 5.8 Coherence of Left Circular Polarization Monterey Bay/La Mesa 



VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

In view of the limited data base and the preliminary 
state of the analysis, only tentative conclusions are 
possible at this point. 

(1) In the frequency range .05 to .4 Hz, the ocean and 
land data exhibit a high degree of coherence. At higher 
frequencies the coherence rapidly diminishes. 

(2) Both polarized and unpolarized signals were observed 
with a preponderance of high degrees of polarization occurring 
at the lower frequencies. The nature of the polarization 
appears to be a function of both time and frequency. 

B. RECOMMENDATIONS 

The first and obvious recommendation is to collect more 
data. This is already in the planning stage. The second and 
also obvious recommendation is to conduct a more thorough ana- 
lysis of the data on hand. This also is already being planned. 
More concrete recommendations concern data processing. 

One 90 -minute real time data tape requires 12 computer 
cylinders of storage space. Thus, for the approximately 
eight hours of data computed and analyzed here, 72 cylinders 
of storage were used. This is too much and is mostly a 
matter of the computer code employed. A little work can 
provide a program that would require less computer storage 



59 



space. This will be a real necessity as a data base of 
weeks or months is compiled for analysis. It is recommended 
that the present code by "streamlined" to use the IBM 
3033 as a calculation "scratch pad." Then a few lines of 
"Computer Job Control Language" can be employed to store 
the calculations on a high-density tape for future use at 
the Naval Postgraduate School or anywhere in the world. 
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APPENDIX A 
TEST PROGRAM 

//JOEM1MB JOB (1592,01 29) ,» FISHER SMC. 1 399 ' , C LASS = C 
//* MAIN OR G=NPGVM1. 159 2P , L INES= ( 60) 

//^FORMAT PR,DDNAME=PL OT.SYSVECTR, DEST= LOCAL 

// EXEC FRTXCLG?,PARM. LKED=’LIST, MAP, XREF' , REGION. GO=2700K 

//SYSUT1 DD UNII=SY SDA, SPACE= (CYL, (8, 8) ) 

//SYSUT2 DD UNIT=SYSDA, SPACE= (CYL, (3,3) ) 

//SYSLIN DD SPACE= (6080 , (80,30) ) ,UNIT=SYSDA 

//FORT. S YS I N DD DSN=MS S. SY S3 . NON IMSL . SOURCE ( FOURT) , DISP= SHR 

// DD * 

C A LARGE NUMBER OF ARRAYS ARE UTILIZED 

C WITH THE INTENTION OF MAINTAINING EACH 

C BIT OF ORIGINAL DATA OR CALCULATION, 

C BEADY FOR RECALL AT ANY TIME FOR FURTHER 

C COMPUTAI TON . EQUIVALENT ARRAYS ARE EMPLOYED 

C EECAUSE THE PLOTTING ROUTINE 'DRAWP' IS 
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c 



NOT ABLE TO HANDLE COMPLEX ARRAYS 



C 



ARRAYS ' IPB ' , ’RTB' , ' ALAB' AND ’TITLE’ ARE 



C 



USED TO GENERATE THE VERSATEC PLOTTER OUTPUT 



C 



OCEAN ARRAYS ARE INDICATED VIA THE SUFFIX 



C 



'O’ AND LAND ARRAYS ARE INDICATED WITH ’L' 



COMP LEX* 8 XXO(3 192) ,YYO(3192) , 

C CSO (8192) ,COO (8 192) , ZXO (3192) , 

C CRO (8192) ,CTO(3 192) ,CUD (3192) , 

C C VO (8192) ,CWO (3 19 2) ,CPO (319 2) , 

C ZYO (819 2) , CO MOL (819 2) ,C0L (8192) , 

C CUL (819 2) , CVL(3 19 2) , CWL (3192) , 

C COPO L(8 1 92) , XXL (819 2) , YYL (819 2) ,ZXL(819 2) f 
C ZYL (8192) ,CSL (8 192) ,BP0 (3 192) , 

C BMO (819 2) , 3PL(S 19 2) ,3ML (819 2) , 

C SPOL (8192) , SM3L (819 2) ,SP0 (819 2) , 

C SPL (8192) , SMO(8 192) ,CTL (3 192) , 
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C CPL (8192) ,CRL(3 192) ,SML(9192) 



DIMENSION TIME (8192) r FRE3 (819 2) , 

C WORK (16384) , FRQ2 (81 92) 

INTEGERS IT B ( 1 2) /I 2*0/ 

REAL *4 RTB(28) /28*0. 0/, 

C CROO (16 3 84) , CFO 0 (16 384) ,2 000 (16 384) , 

C COOO (16 3 84) , CTO 0 ( 16 3 84) , 3VQ0 (16384) , 

C CWOO (16384) ,0111(16384) , 

C CPLL (16 3 84) , CRL L( 16 3 84) , TOLL (16384) , 

C CULL (16384) ,CVL L ( 16 3 84) , 

C CWLL (16 3 84) , COM OLL ( 1 6 384 ) , COPOLL (16384) , 
C SPOO (16384) , S?L L(16384| , 

C SMOO (16384) , SMLL(16384) 

REAL ALAB (22) /' PDXO' , ' PDY D » , » ST 10 ■ , 

C ' ST20» , • ST30 * , • STOO* , 'PHAO 1 , 

C »GAMC» , » PDXL * , ' PDYL 1 , ’ S T 1 L 1 , ' ST 2L ' , 
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C ' ST3L' , ' STOL' , • PHAL 1 , '3AML' , 



C ’ SPL O' , ' SPLL', ' SHOO* , 'SHLL* , ' CPOL • , 1 C MOL ' / 

REAL *8 TITLE (12) 

EQUI VALENCE (TIT LE (1 ) , RTB (5) ) , (CROO(1) , 

C CEO ( 1) ) , (CP30 ( 1 ) , CPO (1) ) , 

C (COOO(1) ,COO (1) ) , (CTOO ( 1| , CTO ( 1) ) , (CUOO (1) ,CUO (1) ) , 

C (CWOO(1) ,CWO (1) ) , (CVOO(1) , CVO (1) ) , (CTLL(1) ,CTL(1) ) , 

C (CPL L ( 1 ) , CPL ( 1 1 ) , (CELL ( 1 1 ,CRL (1) ) , (COLL (1) ,COL (1) ) , 

C (COLL ( 1 ) ,CUL (1) ) , (CVLL(1| ,CVL ( 1) ) , (CWLL (1) ,CWL (1) ) , 

C (COPOLL (1) ,COP3 L ( 1) ) , (COMOLL (1) , 

C COMOL(1 ) ) , (SPOO (1) , SPO (1) ) , 

C (SPLL(1) , SPL (1) ) , (S MOO (1| ,SMO (1) ) , (SMLL (1) ,SML (1) ) 
DATA XXO ,YYO,CSO,COO,ZXO,ZYO,CRO,CTO,CUO,CVO, CWO f 
C XXL, YYL,ZXL,ZYL , CSL, BPO f B MO , BPL , 

C BML, SPOL,SMOL,SPO,SPL,SMO, 

C SML, CTL,CRL,C3L ,CVL f COPOL,CWL, 
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COPO I, CO MOL/278 528* (0.0, 0.0)/ 



DATA TIME, FREQ, FRQ2/2457S * 0 . 0/ 

IFR A ME=8 1 92 
NR=3 0 

FNR = FLOAT (NR) 

P=0 . 0 

DO 7 0 L 1 = 1 , N R 

THE DO LOOP ENDING WITH STATEMENT 70 ENABLES 
THE PROGRAM TO PROCESS A LARGE AMOONT OF DATA 
BY REPEATING THE PROCESS IN BLOCKS. THE DATA 
POINTS FROM EACH RUN THROUGH THE DO LOOP ARE 
ADDED TOGETHER AND EVENTUALLY AVERAGED BY THE 
NUMBER OF RUNS THROUGH THE DO LOOP. 'NR' 
REPRESENTS THE NUMBER OF DATA SEQUENCES TO 
BE AVERAGED. 1 SEQUENCE CURRENTLY EQUALS 



8192 DATA POINTS FOR EACH CHANNEL OR 256 



SECONDS OF DATA. FOR THIS TEST SINE AND 



COSINE WAVES OF 4 HZ NIL BE USED. 

DO 6 0 J J= 1 , I FR A ME 

XXO ( JJ) =SIN (P*. 78 53 98163) 

YYO ( JJ) =COS(P*. 7853 9816 3) 

XXL ( JJ) = SIN (P*. 7853 9815 3) 

YYL ( JJ) =COS (P*. 78 53 98163) 

P=P+ 1 
CONT INUE 

THE FOLLOWING SECTION GENERATES THE TIME 
AND FREQUENCY ARRAYS AND NORMALIZES 
THE INPUT PCM DATA TO VOLTAGE FORM IN 
PREPARATION FOR THE FAST FOURIER TRANSFORM 
TO THE FREQUENCY DOMAIN. 

N=8 1 92 



FN= FLOAT (N) 



DELT AT= 1./32. 



T=FN *DELTAT 
DELT AF=1 ./T 
DO 2 C J = 1 , N 

TIME (J) = DELT AT* FLOAT (J| 

FREQ (J) =DELTAF* FLOAT (J| 

FRQ2 (J) = ALOG 1 0 (FREQ (J) ) 

20 CONTINUE 

C THE NEXT POOR STATEMENTS 

C AN FFT ON THE INPUT TIME 

C DATA. 

CALL FO U RT ( X XO , N , 1 , -1,0,30RK) 
CALL FOURT (YYO, N, 1, -1 ,0, »ORK) 
CALL FOURT(XXL, N, 1, -1 ,3, WORK) 
CALL FOURT (YYL, N, 1, -1,0,ffORK) 
DO 40 K4 = 1 , N 



PERFORM 

SERIES 
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XXO ( K4) =XXO ( K4) /FN 



YYO ( K4) =YYO (K4) /FN 
XXL ( K4) =XXL (K4) /FN 
YYL ( K4) =YYL(K4) /FN 
40 CONTINUE 

C THE NEXT LOOP CALCULATES THE 

C UNNORMALIZED SPECTRAL DENSITIES 

C FOR SINS LE SITE ORTHOGONAL 

C COMPONENTS OF THE GEOMAGNETIC FIELD, 

C THE INDIVIDUAL SITS OROS SPECTRA 

C EETWEEN COMPONENTS, THE INDIVIDUAL 

C ORTHOGONAL COMPONENTS OF RIGHT AND 

C " LEFT CIRCULARLY POLARIZED SPECTRA, 

C AND TWO SITE CROSS SPECTRA FOR RIGHT 

C AND LEFT CIRCULAR POLARIZATION. 

DO 3 C 11 = 1, N 
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ZXO (II) =ZXO(II) + (XXO (II) *C3NJG (XXO (II) ) ) 



ZYO (II) = ZYO (II) + (YYO (II) * CO NJG (YYO (II) ) ) 
CSO ( II) =CSO (II) + (XXO (II) * CO NJG (YYO (II) ) ) 
ZXL (II) =ZXL (II) + (XXL (II) *CONJG (XXL (II) ) ) 

ZYL (II) = ZYL (II) + (YYL (II) *CONJG (YYL (II) ) ) 

CSL (II) =CSL (II) + (XXL (II) * CO N JG (YYL (II) ) ) 

BPO (II) = XXO (II) +((0. ,1.) * YYO (II) ) 

BMO (II) =XXO(II) -( (0. ,1.) *YYO (II) ) 

B?L (II) = XXL (II) + ((0.,1.) * YYL (II) ) 

BML (II) = XXL (II) -( (0. ,1.) *YYL(II) ) 

SPOL (II) =SPOL (I I) + (BPO (II) * CON JG (BPL (II) ) ) 
SMOL (II) = SMOL (I I) + (BMO (II) *CONJG (BML (II) ) ) 
SPO (II) =SPO(II) + (BPO (II) *CONJG (BPO (II) ) ) 
SPL (II) = SPL (II) + (BPL (II) * CO N JG (BPL (II) ) ) 
SMO (II) =SMO (II) + (BMO (II) *CONJG (BMO (II) ) ) 
SML (II) = SML (II) + (BM L (II) *CONJG (BML (II) ) ) 
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3 0 CONTINUE 
70 CONTINUE 



c 


THE NEXT LOOP NORMALIZES THE 


A3 0 V 


c 


SPECTRA AND CALCULATES POWER 


SPECT 


c 


LENS ITIES. 






DO 33 13=1, N 






ZXO (13) =ZXO (13) *T/FNR 





ZYO (13) =ZY0(I3) *T/F NR 
CSO ( 13) =CSO (13) *T/FNR 
ZXL (13) =ZXL(I3) *T/FNR 
ZYL (13) =ZYL(I3) *T/FNR 
CSL (13) =CSL (13) *T/FNR 
SPOL (13) =SPO L (I 3) *T/FNR 
SPO (13) = SPO (13) *T/F NR 
SPL (13) =SPL (13) *T/FNR 
SMO (13) =SMO (13) *T/FNR 
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SML (13) = SML ( 13) *T/FNR 



SHOL (13) = SM0 L (I 3) *T/FNR 
33 CONTINUE 

C THE NEXT LOOP CALCULATES STOKES 

C PARAMETERS 0 THROUGH 3 OF THE 

C INDIVIDUAL SITE ORTHOGONAL 

C COMPONENTS AND THE COHERENCE OF 

C THE PLANAR CIRCULAR POLARIZATION 

C PARAMETERS BETWEEN SITES. 

DO 4 4 14=1, N 

CTO (14) = (ZXO (14 ) +ZYO (14) | *2./T 

CPO ( 14) = (ZXO (14 ) -ZYO (14) ) * 2 ./ (T*CTO (14) ) 

CRO (14) = (4*CS0 ( 14) ) / (CTO (14) *T) 

COO (14) =CSO (14 ) /(CSQRT (ZXO (14) ) *CSQRT (ZYO (14) ) ) 
CWO (14) = AT AN 2 (AIM AG (COO (I 4 ) ) , RE AL (COO ( 14) ) ) 

COO (14) =CSQRT (C 00 (14) *CON JG (COO (14) ) ) 
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CTO (14) =4.34 294 48*CLOG (CTO (14) ) 



CUO (14) =4.34 294 48*CLOG(ZXO (14) ) 

C VO ( 14) =4.34 294 48*CLOG (ZYO (14) ) 

CTL (14) = (ZXL (14 ) +ZYL (14) ) *2./T 

CPL (14) = (ZXL (14) -ZYL (14) ) *2./(T*CTL (14) ) 

CRL (14) = (4*CSL ( 14) ) /(CTL(I4) *T) 

COL (14) =CSL (14) / (CSQRT (ZXL (14) ) *CSQRI (ZYL (14) ) ) 

CTL ( 14) =4. 34 294 48*CLOG(CTL (14) ) 

CUL ( 14) =4.34 294 48*CLOG (ZXL (14) ) 

CVL ( 14) =4.34 294 4 8*CLOG (ZYL (14) ) 

CWL (14) = AT AN 2 ( A IM AG (COL (I 4) ) , REAL (COL (14) ) ) 

COL (14) =CSQRT (COL (14) *CONJG (COL (14) ) ) 

COPO I (14 ) =SPOL ( 14) / (CSQRT (SPO (14) ) *CSQRT (SPL (14) ) ) 
COPO L (I 4 ) =CSQRT (COPOL (14) *CONJG (COPOL (14) ) ) 

COHO 1(14) =SHOL( 14) /(CSQRT (SMO (14) ) *CSQRT (SHL (14) ) ) 
COMO L (14) =CSQRT (COHOL (14) *CONJG (COM3L (14) ) ) 
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SPO (14) =4.34 294 48*CLOG (SPO (14) ) 



SPL (14) = 4. 34 29 448 *C LOG (SPL (14) ) 

S MO (14) =4.34 294 48*CLOG(SMO (14) ) 

SML (14) =4.34 294 48*0 LOG (SML (14) ) 

44 CONTINUE 

C VER3ATEC PLOT OF CALCULATED QUANTITIES 

N PT S =1 0 . /DELTA F + 1 . 

C <N?TS' DETERMINES NUMBER OF POINTS 

C NECESSARY TO PLOT THE 0 TO 10 HERTZ RANGE. 

ITB (2) =0 
IT3 ( 3) =2 0 
ITB (4) = 10 
I TB ( 6) = 1 
ITB (7) = 1 
ITB ( 12) =1 
RTB ( 1) =0 
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RTB ( 2) =0 



RT3 ( 3) = ALAB ( 1) 

READ (5,3000) TITLE 

CALL DRAWP(N?TS , FRQ2,OJOO r ITB, RTB) 

ITB (6) =1 

RTB (3) =ALAB (2) 

READ (5, 3000) TITLE 

CALL DRAW? (NPTS ,FRQ2 ,CV03 , ITB, RTB) 

ITB (6) =1 

RTB (3) = ALAB (3) 

READ (5,3000) TITLE 

CALL DRAWP (NPTS ,FREQ,CPOD, ITB, RTB) 

I TB ( 6) = 1 

RTB ( 3) = A LAB (4) 

READ (5, 3000) TITLE 

CALL DRAWP (NPTS , FREQ , CROD , ITB , RTB) 
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ITB (6) = 1 



RT3 ( 3) = A LAB (5) 

READ (5, 3000) TITLE 

CALL DRAWP(NPTS , FRE Q, CR03 ( 2 ) , IT E, RT B) 

ITB (6) = 1 

RTB ( 3) = AL AB ( 6) 

READ (5,3 000) TITLE 

CALL DRA WP (N PIS , FRQ 2 , CTOO , I TB, RTB) 

ITB ( 6) = 1 

RTB ( 3) = A LAB ( 7) 

READ (5, 3000) TITLE 

CALL DRAW? (NETS , FREQ, CWOD , IT3, RTB) 

IT3 (6) = 1 
RTB (3) = ALA3 (8) 

READ (5,3 000) TITLE 

CALL DRAWP(N?TS , FREQ, COOD , ITB, RTB) 
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I TB ( 6) = 1 



RTB (3) =ALAB(9> 

READ (5 , 3 000) TITLE 

CALL DRAWP (N PI'S , FRQ2, CULL , ITB, RTB) 

ITB (6) =1 

RTB (3) =A LAB (10) 

READ (5, 3000) TITLE 

CALL DRAWP(NPTS,FRQ2,CVLL,ITB,RTB) 

ITB (6) = 1 

RTB ( 3) =ALAB (11) 

READ (5, 3 0 00) TITLE 

CALL DRAWP (NPT3 , FREQ, CPLL , ITB, RTB) 

ITB (6) =1 

RTB ( 3) =A LAB ( 12) 

READ (5, 3000) TITLE 

CALL DRAWP (NPTS , FRE Q , CRLL , I TB, RTB) 
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ITB (6) = 1 



RTB ( 3) =A LAB (13) 

READ (5, 3000) TITLE 

CALL DRA WP (NPTS , FREQ, CELL (2) , ITB, RI3) 

ITB (6) = 1 

RTB (3)=ALAB(14) 

READ (5, 3000) TITLE 

CALL DRA WP (NPTS , FRQ2 , CTLL , ITB, RTB) 

ITB (6) =1 

RTB ( 3) = A LAB ( 15) 

READ (5, 3000) TITLE 

CALL DRA WP (NPTS , FRE Q , CWLL , ITB , RT 3) 

ITB (6) = 1 

RTB (3)=ALAB(16) 

READ (5, 3 000) TITLE 

CALL DRA WP (NPTS , FREQ, COLL , ITB, RTB) 
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ITB (6) =1 



RTB (3) =ALAB( 17) 

READ (5,3000) TITLE 

CALL DR AWP (NPIS , FRQ 2 , SPOO , I TB , RTB) 

ITB (6) =1 

RTB (3) =ALAB (18) 

READ (5, 3000) TITLE 

CALL DRA WP (NPIS , FRQ2 , SPLL , I TB, RTB) 

ITB (6) = 1 

RTB (3) =ALAB( 19) 

READ (5, 3000) TITLE 

CALL DR AWP (NPIS , FRQ 2 , SMOD , ITB, RTB) 

ITB (6) =1 

RTB ( 3) = A LAB (20) 

READ (5, 3000) TITLE 

CALL DRA WP (N PTS , FRQ2 , SHLL , I TB , RTB) 
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ITB (6) = 1 



RT3 (3) =ALAB (21) 

READ (5, 3000) TITLE 

CALL DR A WP (NPTS ,FREQ, COPOLL , IT E, RTB) 

ITB (6) =1 

RTB ( 3) = ALAB (22) 

READ (5, 3000) TITLE 

CALL DR A WP (NPTS , FREQ, COMOLL , ITB, RTB) 
3000 FORMAT (6A8) 

STOP 

END 

/* 

//GO. SYS IN ED * 

TEST POWER SPECTRAL DENSITY OF SIN 4HZ 
IN XXO,30 A VGS , (REF DB YS LOG FREQ) 

TEST POWER SPECTRAL DENSITY OF COS 4HZ 
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IN YYO,30 A VGS , (REF DB VS LOG FREQ) 
TEST STOKES 1/STOKES 0 OF SIN 4HZ 
IN XXO AND COS 4HZ IN YYO, 30 AVGS 
TEST STOKES 2/STOKES 3 OF SIN 4HZ 
IN XXO AND COS 4HZ IN YYO, 30 AVGS 
TEST STOKES 3/STOKES 0 OF SIN 4HZ 
IN XXO AND COS 4HZ IN YYO, 30 AVGS 
TEST STOKES 0 OF SIN 4 HZ IN XXO AND 
COS 4HZ IN YYO, 30 AVGS, (REF DB VS LOG) 
TEST PHASE OF SIN 4 HZ IN XXO 
AND COS 4HZ IN YYO, 33 AVGS 
TEST COHER OF SIN 4 HZ IN XXO 
AND COS 4HZ IN YYO, 33 AVGS 
TEST POWER SPECTRAL DENSITY OF SIN 4HZ 
IN XXL, 30 AVGS, (REF DB VS LOG FREQ) 
TEST POWER SPECTRAL DENSITY OF COS 4HZ 
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IN YYL, 30 AVGS. (REF D3 VS LOG FREQ) 

TEST STOKES 1/STOKES 0 OF SIN 4HZ 
IN XXL AND COS 4HZ IN YYL, 30 AVGS 
TEST STOKES 2/STOKES 0 OF SIN 4HZ 
IN XXL AND COS 4HZ IN YYL, 30 AVGS 
TEST STOKES 3/STOKES 0 OF SIN 4HZ 
IN XXL AND COS 4HZ IN YYL, 30 AVGS 
TEST STOKES 0 OF SIN 4 HZ IN XXL AND 
COS 4HZ IN YYL, 30 AVGS, (REF DB VS LOG) 

TEST PHASE OF SIN 4 HZ IN XXL 
AND COS 4HZ IN YYL, 30 AVGS 
TEST COHER CF SIN 4 HZ IN XXL 
AND COS 4HZ IN YYL, 30 AVGS 

TEST RT CIRC POLARIZATION PSD OF SIN 4 HZ IN 
XXO AND COS 4HZ IN YYD , (DB VS LOG), 30 AVGS 
TEST RT CIRC POLARIZATION PSD OF SIN 4HZ IN 
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XXL AND COS 4HZ IN YYL , (DB VS LOG) ,30 AVGS 



TEST LEFT CIRC POLARIZATION PSD OF SIN 4HZ 
IN XXO AND COS 4HZ IN YYO (DB VS LOG) ,30 AVGS 
TEST LEFT CIRC POLARIZATION PSD OF SIN 4HZ 



IN XXL AND 


COS 4HZ IN YYL (DB VS LOG), 30 AVGS 


TEST COHER 


FT CIRC POLARIZATION OF SIN 4HZ 


IN XXO, XXL 


AND COS 4HZ IN YYO, YYL, 30 AVGS 


TEST COHER 


LEFT CIRC POLARIZATION OF SIN 4 HZ 



IN XXO, XXL AND COS 4HZ IN YYO, YYL, 30 AVGS 
/* 

// 
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Figure A.l Test Ocean x Coil Power Spectral Density 
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Figure A. 2 Test Ocean y Coil Power Spectral Density 
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Figure A. 3 Test Ocean Stokes Parameter Zero 
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Figure A. 4 Test Ocean Stokes Parameter One 
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Iigure A. 5 Test Ocean Stokes Parameter Two 
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Figure A. 6 Test Ocean Stokes Parameter Three 
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FREQUENCY (HZ) 
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FREQUENCY (HZ) 
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Figure A. 9 Test Ocean Right Circular Polarization Power Spectral Density 
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Figure A. 10 Test Ocean Left Circular Polarization Power Spectral Density 
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Figure A. 11 Test Land x Coil Power Spectral Density 




Kl 




>< 

cj 

w 

:=> 

ex 

w 

e* 

tlM 

CJ 

o 

^4 



(ZH/ZvvlN d3U) sa 



94 



Figure A. 12 Test Land y Coil Power Spectral Density 
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Figure A. 13 Test Land Stokes Parameter Zero 
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Figure A. 14 Test Land Stokes Parameter Two 
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I'igure A. 15 Test Land Stokes Parameter Three 
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Figure A. 16 Test Land Phase of x and y Coils 
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FREQUENCY (HZ) 

Figure A. 17 Test Land Coherence of x and y Coils 
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Figure A. 18 Test Land Right Polarization Power Spectral Density 
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Figure A. 19 Test Land Left Polarization Power Spectral Density 
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Figure A. 20 Test Ocean Coherence of Right Circular Polarization 
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Figure A. 21 Test Ocean Coherence of Left Circular Polarization 




APPENDIX B 

MASS STO RAGE PROGRAM 

//JOEB1MB JCB ( 1592,01 29) ,» FISHER SMC. 1399' ,CLASS = F 
//♦MAIN OR G =N?G VMl . 1592P 
// EXEC FORTXCLG 

//FORT . S YS I N DD DSN=HS S. SYS3. NONI MSL. SOURCE ( FOURT) , DISP= S HR 
// DD * 

C THE ARRAY 'IN’ WILL BE USED TO 

C RECEIVE THE DATA PASSED FROM 

C THE SUBROUTINE ‘READ’ AND THEN 

C TRANSFERRED TO THE APPROPRIATE 

C XXX OR YYY ARRAY. 

IN TEG ER*2 IN(1S) 

COMPL EX* 8 XXX (81 92) , YYY (31 92) 

DIMENSION FREQ (9 192) , TIME (8 192) , WORK (16384) 

DATA XXX, YYY/153 84*( 0.0,0. 0) / 

DATA TIME, FREQ/1 6384*0.0/ 
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c 



THE FOLLOWING SECTION READS THE 



C FIRST SIX SECONDS OF COMPUTER TAPE 

C AND DISCARDS THIS DATA. 

ITL= 1 92 

DO 55 JJ=1 r ITL 

CALL RD (20, IN, 20 0, IR EC, IRR) 

55 CONTINUE 
IFRAME=8192 
NR =2 0 

FNR= FLOAT (NR) 

DO 7 0 L 1 = 1 , NR 

C THE NEXT LOOP READS THE COMPUTER 

C TAPE USING THE PROGRAM PROVIDED 

C BY MR. TIM STANTON OF NAVAL POST- 

C GRADUATE SCHOOL. 

DO 6 0 J J= 1 , IFRAM E 
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CALL ED(20,IN,13 00,IREC,IRR) 



XXX(JJ) = IN (2) 

YYY( J J) =IN (3) 

CONTINUE 
N=81 9 2 
FN=FL OAT (N) 

DELTA 1= 1 ./32. 

T=FN* DELT AT 
DELTA F=1 . /T 
DO 20 J=1 ,N 

TIME ( J) =DELT AT*F LOAT ( J) 

FREQ (J) = D ELTAF*F LOAT (J) 

THE NEXT 4 STEPS CONVERT THE 

DATA TO VOLTAGE AND ENSURES THAT 

NO ERRONEOUS DATA HAS BEEN INTRODUCED 



INTO THE ARRAYS 



XXX (J ) = (XXX(J) -2 045. 5) *1 0. /2045. 5 



XXX ( J) =RE AL (XXX ( J) ) 

YYY (J) = (YYY(J) -2 045. 5 ) * 1 0. /2045 . 5 
YYY ( J)=RE AL (YYY ( J) ) 

20 CONTINUE 

CALL FOURT (XXX, N ,1,-1 ,0, WORK) 

CALL FOURT (YYY, N , 1 ,-1 , 0, WORK) 

DO 40 K4=1,N 
XXX (K 4) =XXX ( K 4> / FN 
YYY (K4) =YYY(K4) /FN 
40 CONTINUE 

C THE NEXT LOOP APPLIES THE SYSTEM 

C TRANSFER FUNCTION TO THE TRANSFORMED 

C FREQUENCY DOMAIN DATA. THE TRANSFER 

C FUNCTION CONVERTS VOLTS TO NANOTESLAS. 

DO 9 L=1 , N 
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FRQ= F BEQ (1) 



1 IF (FRQ.LE.15.) 30 TO 2 

XXX ( L ) =XXX (L) / ( 1 05-5-3. 1 4 * FRQ) 

YYY (L) =YYY (L) / ( 1 3 1 . 3 2- 7 . 59 8*FRQ) 

GO TO 8 

2 IF (FHQ. LS. 10. ) GD TO 3 

XXX ( L ) =XXX (L) / (5 . 958 *FRQ-3 0 . 9 7) 

YYY ( L ) =Y Y Y (L) / ( 7 . 166*FR2-39. 99) 

GO TO 3 

3 IF (FRQ.LE.7.5) GO TO 4 

XXX {L ) =XXX (L) / ( 3 . 492*?R2“5 • 31) 

YYY ( L ) =Y Y Y (L) /(4 .252*FRQ-1 0. 85) 

GO TO 8 

4 IF (FRQ.LE. 5.) GO TO 5 

XXX ( L ) = X X X (L) / (2 .6311 *FRQ*0. 14667) 
YYY (L)=YYY (L) /(3 . 0 1 2* FRQ- 1 . 5 5 ) 
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GO TO 8 



5 IF (FRQ. LE.3. ) 30 TO 6 

XXX ( L) =XXX (L) / (2 .631 1*FRQ*0. 146 67) 

7 YYY(L)=YYY (L) / ( 2 ,702*FRQ) 

GO TO 8 

6 XXX(L)=XXX (L) /(2 .7 2* FRQ) 

GO TO 7 

8 CONTINUE 

9 CONTINUE 

C THE NEXT WRITE STATEMENTS SEND 

C THE CONVERTED DATA TO MSS 

C FOR FUTURE MANIPULATION AND RECALL. 

C THIS NEXT SET OF STATEMENTS ALSO 

C PROVIDES THE USER A DIAGNOSTIC. 

WRITE (21) XXX 

WRITE (6,*) XXX (1) ,XXX (819 2) 
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WRITE (21) YYY 



WRITE <6,*)YYY(1) , YYY (8192* 

70 CONTINUE 
ENDFILE 21 
STOP 
END 

SUBROUTINE P D (I UN , 10 , IRS, IREC, IRQ) 



C 



C 



THIS PROCEDURE FURNISHED BY MR. TIM STANTON 



C 



DEPARTMENT OF OCEANOGRAPHY 



C 



c 



READ DATA FROM IUN, ALLIGN , CHECK & RETURN 



C 



C 



IUN =TA PE NUMBER. EG 20 



C 



10= INTEGER* 2 ARRAY, 16 LONG 



C 



(VALUES 0-4095 , SUBTRACT 2048) *5/2028 . GIVES VOLTAGE 
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c 


IR S = NUMBER DF 


PE5INCS ALLOWED (ERRORS) 


c 


IREC= COUNTER 


OF RECORDS (FRAMES OF DATA) 


c 


BLOCK 512 


BITS, 32 BITS = RECORD 


c 


800 BPI TA 


PE UNLABLED 


c 


IRQ = NUMBER OF 


ACTUAL RESINCS (ERRORS) 


c 






c 








INTEGER * 2 ID ( 


16) , IP ( 1 6) 




DATA IRR /0/ 






IF (IREC.EQ.O) 


15=3 




I ER= 0 




20 


FORMAT ( 1 6A2 ) 






IF (IS.NE.O) SD 


TO 50 



READ (IUN,20,END=900) IP 
IREC=IREC+1 
40 IS=IS+1 
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IF (IS. LT. 17) GO TD 50 



READ (IUN,20,END=900) IP 
IS = 1 

IREC=IREC + 1 

50 ICH=IMASK (IP(IS) ,3,0) +1 
C WRITE (6,55) IC H, IS , IUN, IREC 

55 FORMAT (' RESYNCING ICH, IS ,IUN , IREC ’,413) 

C 

IF (ICH.NE.1) GO TO 40 
DO 100 1=1,16 
IO(I)=ISHIFT (IP (IS) ,4) 

ICH= I MASK (IP (IS) ,3,0) +1 
IF (ICH.EQ.I) GO TO 80 
IER= I ER+ 1 

WRITE (6,70) IUN, IREC ,I,ICH,IER 
70 FORMAT ( 1 UNIT* , 13 ,' RECORD ', 16 ,’ CHAN S DATA CH ',214, 
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'ERRORS ',17) 



80 IS=IS+1 

IF (IS.LT.17) 30 TO 100 
READ (ION, 20, EM D=900) IP 
IS=1 

IREC =IREC+1 
100 CONTINUE 
C 

IF ( IER. EQ.9) 30 TO 150 
I RR= IRR + 1 

IF (IRR. LT. IRS) GO TO 120 
WRITE (6,110) 

110 FORMAT (’1 STOPPED IN SUB RD BECAUSE 
C OF IRR.GT.' ,15, ’ AT L110') 

IRQ=IRR 

STOP 
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120 CONTINUE 

WRITE (6,130) IREC, IRR 

130 FORMAT ('RESYNC AT FRAME ' , 16 , ' W IT H TOTAL ERRORS', 17) 
ISR= 0 
IRQ = IRR 
GO TO 50 
150 CONTINUE 
RET U 5N 

900 WRITE (6,910) I UN, I REC 

910 FORMAT ('1 EN D OF UNIT ',13,' AT REC ',17) 

STOP 

END 

FUNCTION ISHIFT (IN,NPLC) 

C RETURNS SHIFTED VALUE OF 1*2 WORD IN 

C -VE LE FT , + VE RIGHT SHIFT 

C 
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INTEGER * 2 IN 



IP= I N 

IF (IP.LT.O) I?=IP* 65536 
IF (NPLC.LT.O) GO TO 30 
I SH I FT=I P/ (2 **I ABS(NPLC) | 

RET U RN 

30 I SH I FT=I P* (2**1 ABS(NPLC) ) 

IF (ISHIFT.GT. 65535) ISHIFT = MOD (ISHIFT , 65536) 

RET U RN 
END 

FUNCTION IMASK (IN, IBL, IBR) 

C MASK 1*2 WORD IN OUTSIDE BITS IBL S IBR 

C 

INTEGER * 2 IN, 10 
10= IN 

IF (IBR.FQ.O) 30 TO 50 
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IT= I SHIFT (IN , I B R) 



10=1 T 



50 IP=ISHIFT (I0,IBL-15-IBR) 



IO=IP 



I MAS K= IS HIFT (ID , 15-IBL) 



RETURN 



END 



/* 



//GO. FT2 1F0 01 
// 

// 

// 

//GO. FT20FO 01 
// 

// 

// 



DD UNIT=3330V,MSV3P=PUB4A,DISP=(NEW,CATLG) , 

D SN=M SS . S1592.GMDI1 A, 

DCB= ( REC FM=VBS,BLISIZE=4096,LRECL=4092) , 

SPACE= (0 YL, (8 ,8) ) 

DD UN 11=3400-4, VOL = SER = GMDT1A,DISP= (OLD, PASS) , 
L A3EL= ( 1 , NL, , IN) , 

DCB= (RECFM=FB,LREGL=32, BLKSIZ E=5 12,DEN = 2) 
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APPENDIX C 

\ MAIN PROGRAM 

//JOEK1MB JOB ( 1592,01 29) , ' FISHER S MC. 1 399 » , CLASS= F 

//♦MAIN OR G=NPGVM1 . 159 2P , L INES= (60) 

//♦FORMAT PR,DDNAME=PLOT.SYSVECTR, DEST= LOCAL 

// EXEC FRTXCLGP,PARM. LKED= ' LIST , M AF , XREF' , R SGI ON . GO = 270 OK 

//SYSUT1 DD UNIT=SI SDA, SPACE= (CYL, (8, 8) ) 

//SYSUT2 DD UNIT=SY SDA, SPACE= (CYL, (8,8) ) 

//SYSLIN DD SPACE= (6080 , (80,30) ) ,UNIT=SYSDA 

//FORT . S YS I N DD DSN=MS S. SYS3 . NON IMS! . SOURCE (FOURT) , DI SP= SHR 

// DD * 

C DETAILED DESCRIPTIONS OF THE 

C VARIOUS STEPS AND LOOPS ARE 

C CONTAINED IN APPENDICES A AND B. 

INTEGER*2 IN (15 ) 

COMP IEX*8 XXO(8 192) , YYO (3 192) , 

C CSO (8192) ,COO(8 192) , ZX0 (3 19 2) , 
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C CRO (819 2) ,CTO (8 19 2) ,CUO (3192), 

C CVO (8192) ,CWO(8 192) ,CPO(3192) , 

C ZYO ( 8192) , CO MOL (819 2) ,COL (8192) , 

C CUL (819 2) ,CVL(8 192) ,CWL (319 2) , 

C COP OL(819 2) , XXX (8192) ,YYY (8192) , 

C XXL ( 8192) ,YYL(8 192) , ZXL (3 192) , 

C ZYL (819 2) ,C5L(8 192) ,BPD (819 2) , 

C BMC (8192) ,BPL(9 192) ,BML (3 19 2) , 

C SPOL (8192), SMOL (8192) ,SPD (8192) , 

C SPL (8192) ,SM0(3 192) ,CTL (8192) , 

C CPL (8192) ,CRL(3 192) ,SML(3192) 

DIMENSION TIME ( 8192) , FRE2 (8 192) , 

C WORK (16384) , FR2 2 (81 92) 

INTEGER*4 ITB ( 1 2) /I 2 *0/ 

REAL *4 RT B (2 8) /2 8*0. 0/, CROO (16 384) , 
C CPOO (16384) , CUD0(16384) , 
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C COOO (16384) , CT 00(16384) ,CVOO (16384) , 

C CWOO (16384) , dL L( 16384} , 

C CPLL (16384) ,CRLL(1S384) , COLL (16384) , 

C CULL (16384) , CVLL (16384) , 

C CWLL (16384) , CONOLL ( 1 6384> , COPOLL ( 16384) , 

C SPOO (16384) , SPLL(15384) , 

C SMOO (16384) , SML L(16384) 

REAL ALAB (22) /' PDXO» , ' PDYO’ , 1 ST 10* , 

C ' ST20» , ' ST30» , ' STOO' , 1 PHAO' , 

C • GAMO* , ’ PDXL ' , ' PDYL * ,»ST1L»,»ST2L , f 
C »ST3L' , ' STOL',' ?HAL» ,'3AML* , 

C » SPL O' , • SPLL ' , • SMOO 1 , ’ SMLL ' , • CPOL* , • CMOL' / 

REA L *8 TITLE (12) 

EQUIVALENCE (TITLE (1 ) , RTB{5) ) , (CR00(1) , 

C CRO ( 1) ) , (CP00(1 ) , CPO (1>) , 

C (C000(1) ,C00 (1) ) , (CTOO ( 1 ) , CTO (1) ) , (CU00(1) f CUO (1) ) , 
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C (CW00(1) ,CWO (1|), (C VOO ( 1 * , CVO ( 1) ) , (CTLL (1) ,CTL(1) ) , 

C (CPL L(1 ) , CPL (1) ) , (CRLL ( 1 1 , CRL(1) ) , (COLL (1) , COL (1) ) , 

C (CUL L(1 ) ,CUL (1) ) , (CVLL ( 1 ) ,CVL(1) ) , (CWLL ( 1 ) , CWL ( 1) ) , 

C (COPCLL ( 1) ,COPO L (1) ) , (COMOLL (1) , 

C COMOL(I) ) , (SPDO (1) ,SPO (1) ) , 

C (SPL L ( 1 ) , SPL (1) ) , (SMOG (1) , SMO (1) ) , (SMLL(1) , SML(1) ) 

DATA XXO , YY 0 , C S 0,C0 0,ZX0, ZYO, CRO, 

C CTO, CUO, CVO, CWO ,XXX ,YYY, 

C XXL, YYL,ZXL,ZYL , CSL , BPO , B MO , B PL , 

C B ML, SPOL , SMOL, S PO , SPL, SMO , 

C SML , CTL , CRL, COL ,CVL , COPOL , CWL , COPOL , 

C COMO L/29 491 2 *(0.0, 0. 0) / 

DATA TIME, FREQ, FRQ2/24575 *0. 0/ 

ITL= 192 

DO 55 J J= 1 , ITL 

CALL RD (20, IN, 2 00,IR2C,IRR) 
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55 CONTINUE 

IFR A ME = 8 1 92 
NR= 2 0 

FNR = FLOAT (NR) 

C THE FIRST STATEMENTS OF LOOP 70 

C RECALLS DATA PREVIOUSLY STORED IN 

C THE IBM 3033 MASS STORAGE SYSTEM. 

DO 7 C L 1 = 1 , NR 
READ (21) XXX 
READ (21) YYY 
DO 43 II=1,IFRAME 
XXL (II) =XXX (II) 

YYL (II) = YYY (II) 

43 CONTINUE 

DO 6 0 JJ=1, IFRAME 

CALL RD (20,IN,1 000, IREC,IRR) 
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XXO ( JJ) =IN (2) 



YYO ( JJ) = IN (3) 

60 CONTINUE 
N=8 1 92 
F N=F LOAT (N) 

DELTAT=1./32. 

T = FN *DELT AT 
DELT AF= 1 ./T 
DO 20 J= 1 , N 

TIME (J) = DELT AT * FLOAT (J) 

FREQ (J) =DELTAF* FLOAT (J) 

FRQ2 (J) = ALOG10 ( FREQ (J) ) 

XXO (J)= (XXO(J) - 2045.5) *5./2045.5 
XXO ( J) =REAL (XXO (J) ) 

YYO ( J) = (YYO (J) - 2045 .5) *5. /2 045. 5 
YYO ( J) =R EAL (YYO (J) ) 
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20 CONTINUE 



CALL FOURT(XXO, N, 1, -1 ,0 , WORK) 

CALL FOURT{YYO, N, 1, -1,0, WORK) 

DO 40 K4 = 1 , N 
XXO ( K4) =XXO(K4) /FN 
YYO (K4) =YYO (K4> /FN 
40 CONTINUE 

DO 9 L= 1 , N 
FRQ = FREQ (L) 

1 IF { F RQ. LE.15.) 30 TD 2 

XXO (L) =XXO (L) /(-81. 36 +17.066*FRQ - 0 . 5389 *FRQ*FR Q) 
YYO (L) =YYO (L) /( -63. 435 +1 1 . 71 06*FRQ-0. 23279*FRQ*F RQ) 
GO TO 8 

2 I F ( F RQ . LE.10.) 3 0 TD 3 

XXO (I) =XXO (L) /( -18. 486 + 4.9029*FRQ) 

YYO (L) =YYO(L) /(-63. 435 +1 1.7106*FRQ -0 . 2 3 27 9*FRQ* FRQ) 
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GO TO 8 



3 IF(FRQ. IE. 7. 5) GO TO 4 

XXO (L) =XXO (L) /( 3. 947*FRQ-9. 368) 

YYO (I) = Y YO (L) /( 4. 295*FRQ- 1 1.837) 

GO TO 8 

4 IF(F HQ. LE. 5. ) GO TO 5 

XXO (L) =XXO (L) /( 2.8105*FFQ-0. 56 57) 

YYO (L) =YYO (L) / ( 2. 80 0*FRQ-. 5 0) 

GO TO 8 

5 XXO (L) =XXO (L) /( 2. 723*FRQ) 

YYO (L)=YYO (L)/(2.71 7*FRQ) 

GO TO 8 

8 CONTINUE 

9 CONTINUE 

DO 3 0 11 = 1, N 

ZXO (II) =ZXO (II) + (XXO (II) * CONJG (XXO (II) ) ) 
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ZYO (II) = Z YO (II) + (YYO (II) *CONJG (YYO (II) ) ) 



CSO(II) = CSO(II) +(XXO(II) * CO NJG (YYO (II) )) 

ZXL (II) =ZXL (II) + (XXL (II) *C0 N JG (XXL (II) ) ) 

ZYL (II) = ZYL (II) + (YYL (II) * CO NJG (YYL (II) ) ) 

CSL (II) = CSL ( II) + (XXL (II) * CO NJG (YYL (II) ) ) 

3 PO (II) =XXO(II) + ((0.,1.) * YYO (IT) ) 

BMO (II) =XXO (II) - ( (0. ,1.) * YYO (II) ) 

BPL (II) =XXL ( II) + ( (0. f 1.) * Y YL (II) ) 

BML (II) = XXL (II) -( (0. ,1.) *YYL (II) ) 

SPOL (II) =SPOL (I I) + (BPO (II) * CON JG (BPL (II) ) ) 
SMOL (II) = SMOL (II) + (3H0 (II) *CONJG(BML (II) ) ) 
SPO (II) =SPO (II) + (BPO (II) *CONJG (BPO (II) ) ) 
SPL (II) = SPL (II) + (BP L (II) * CON JG (BPL(II) ) ) 
SMO (II) = SMO (II) + (BMO (II) *CONJG (BMO (II) ) ) 
SML (II) =SML(II) + (BML (II) *CONJG (BML (II) ) ) 

30 CONTINOE 
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70 CONTINUE 



DO 33 13=1, N 
ZXO (13) =ZXO (13} *T/FNR 
ZYO (13) =ZYO ( 13} *T/FNR 
CSO (13) =CSO (13} *T/FNR 
Z XL (13) =ZXL (13) *T/F NR 
ZYL ( 13) =ZYL (13) *T/FNR 
CSL (13) =CSL (13) *T/FNR 
SPOL (13) =SPOL (I 3) *T/FNR 
SPO ( 13) =SPO (13) *T/FNR 
SPL (13) =SPL (13) *T/FNR 
SMO (13) = SMO (13) *T/FNR 
SML (13) =SML (13) *T/FNR 
SMOL (13) = SMOL (I 3) *T/FNR 
3 3 CONTINUE 

DO 44 14=1, N 
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CTO (14) = (ZXO (14 ) +ZYO (14) | *2. /T 



CPO (14) = (ZXO (14 ) -ZYO (14) ) * 2 . / <T*CTO (14) ) 

CRO ( 14) = (4*CS0 ( 14) ) / (CTO (14) *T) 

COO (14) =CSO (14) /(CSQRT (ZXO (14) ) *CSQRT (ZYO (14) ) ) 
CWO (14) = ATAN2 (A IKAG (COO (14) ) , RE AL (COO (14) ) ) 

COO ( 14) =CSQRT (COO (14) *C05J JG (COO (14) ) ) 

CTO (14) =4.34294 48*CLOG (CTO (14) ) 

CUO (14) =4.34294 48*CLOG (ZXO (14) ) 

CVO ( 14) =4. 34 29 4 48*CLOG (ZYO (14) ) 

CTL (14) = (ZXL (14) +ZYL (14) ) *2./T 

CPL (14) = (ZXL (14 ) -ZYL (14) ) * 2 . / (T*CTL (14) ) 

CRL ( 14) = (4*CSL ( 14) ) / (CTL (14) *T) 

COL (14) =CSL (14) /(CSQRT (ZXL (14) ) *CSQRT (ZYL (14) ) ) 
CTL (14) =4.34 294 48*CLOG (CTL (14) ) 

CUL ( 14) =4.34 294 48*CLOG(ZXL (14) ) 

CVL (14) =4.34 294 48*CLOG(ZYL (14) ) 
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CWL (14) = ATAN2 (A IMAG (COL (14) ) , RE AL (COL (14) ) ) 



COL (14) =CSQRT (COL (14) *C0NJG (COL (14) ) ) 

COPO 1(14) =S?OL ( 14) / (CSQRT (SPO (14) ) *CSQRT (SPL(I4) ) ) 
COPOL(I4) =CSQRT (COPOL (14) *CONJG (COPOL (14) ) ) 

COMO L (I 4) =SMOL (14)/ (CSORT (SMO (14) ) *CSQRT (SML (14) ) ) 
COMO L (1 4 ) =CSQRT (COMOL (14) *CONJG (COMOL (14) ) ) 

SPO ( 14) =4.34 29 4 48*CLOG(SP3 (14) ) 

S PL (14) =4.34 294 48*CLOG ( SPL (14) ) 

SMO (14) =4.34 294 4 8*C LOG (SMO (14) ) 

SML ( 14) =4.34 29 4 48*C LOG (SML (14) ) 

44 CONTINUE 

NPTS =10 . /DELTA F + 1 . 

ITB ( 2) =0 
ITS ( 3) =20 
I TB ( 4) = 1 0 
ITB (6) = 1 
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ITB (7) =1 



ITB ( 12) = 1 
RTB ( 1) =0 
RTB ( 2) =0 
RTB (3) =ALAB (1) 

READ (5, 3000) TITLE 

CALL DRA WP (NPTS , FRQ2 ,CUOD , ITB, RTB) 

ITB (6) =1 

RTB (3) =ALAB(2) 

READ (5 r 3000) TITLE 

CALL DRA WP (NPTS , FRQ2 , CVOD , ITB, RTB) 

ITB (6) =1 

RTB (3) =ALAB (3) 

READ (5 r 3 0 00) TITLE 

CALL DRA WP (NPTS ,FREQ,CPOD f IT3,RTB) 
ITB (6) = 1 
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RTB ( 3) = A LAB (4) 



READ (5, 3000) TITLE 

CALL DRAWP(NPTS , FREQ, CR OD , I TB, RTB) 

ITB (6) =1 
RTB (3) =ALAB(5) 

READ (5, 3 000) TITLE 

CALL DRAWP (NPTS , FREQ, CROO (2) , ITB, RTB) 

ITB (6) =1 

RTB ( 3) = A LAB (6) 

READ (5, 3000) TITLE 

CALL DRAWP (NPTS , FRQ2 , CTOD , IT 3 , RTB) 

ITB (6) =1 
RTB (3) =ALAB (7) 

READ (5,3000) TITLE 

CALL DRAWP(NPTS , FREQ, C WOO , ITB, RTB) 

ITB (6) =1 
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RT3 ( 3) = ALAB (8) 



READ (5, 3000) TITLE 

CALL DRA WP (NPTS ,FREQ,COOO,ITB,RTB) 

ITB (6) =1 

RTB (3) = A LAB (9) 

READ (5, 3000) TITLE 

CALL DR A WP (NPTS , FRQ2 , CULL , I TB , RTB) 

ITB (6) = 1 

RTB (3) =ALAB(10) 

READ (5, 3000) TITLE 

CALL DRA WP (NPTS , FRQ2 , CVLL , IT3, RTB) 

ITB (6) =1 

RTB (3) =ALAB (11) 

READ (5, 3000) TITLE 

CALL DRA WP (NPTS , FREQ, CPLL, ITB, RTB) 
ITB (6) = 1 
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RTB (3) =A LAB (12) 



READ (5 # 3000) TITLE 

CALL DR A WP ( NPT S , FREQ, CRLL , ITB, RTB) 

ITB (6) =1 
RTB(3)=ALA3(13) 

READ (5, 3000) TITLE 

CALL DRAWP(NPTS , FREQ, CRL L (2) , ITB, RTB) 

ITB (6) =1 

RTB ( 3) = A LAB (14) 

READ (5, 3000) TITLE 

CALL DRAWP(HPIS,FRQ2,CTLL,ITB,RTB) 

ITB (6) = 1 

RTB ( 3) = AL AB ( 1 5| 

READ (5, 3000) TITLE 

CALL DRAWP (NPTS , FREQ , CWL L , I TB , RTB) 

I TB ( 6) = 1 
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RTB (3) = A L AB (16) 



READ (5, 3000) TITLE 

CALL DRA WP (NPTS , FREQ, COLL, ITB, RTB) 

ITB (6) =1 

RTB (3) =ALAB (17) 

READ (5, 3000) TITLE 

CALL DRA WP (NPTS ,FRQ2,SPOD,ITB,RTB) 

ITB (6) = 1 

RTB ( 3) = AL AB ( 1 8) 

READ (5 , 3 000) TITLE 

CALL DRAWP(NPTS,FR22,SPLL r ITB,RTB) 

I TB ( 6) = 1 

RTB ( 3) =ALAB (19) 

READ (5, 3000) TITLE 

CALL DRA WP (NPTS , FRQ2 , SMOD , ITB , RTB) 
ITB (6) =1 
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RTB (3) = A LAB (20) 



READ (5,3000) TITLE 

CALL DRAWP(NPrS , FRQ2 , SMLL , ITB, RTB) 

ITB (6) =1 

RTB (3) =ALAB(21) 

READ (5, 3000) TITLE 

CALL DR A WP (NPIS , FREQ, C3P0LL , ITB, RTB) 

ITB (6) =1 

RTB ( 3) =ALAB (22) 

READ (5, 3000) TITLE 

CALL DR A WP (NPTS , FREQ , COMO LL , IT B, RTB) 
3000 FORMAT (6A8) 

STOP 

END 

SUBROUTINE RD(I UN, IO, IRS, IREC,IRQ) 

C 
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c 



THIS PROCEDURE FURNISHED BY HR. TIM STANTON, 

C DEPARTMENT OF OCEANOGRAPHY. 

C 

C READ DATA FROM ION, ALLIGN , CHECK & RETURN 

C 

C ION =T A PE NUMBER, EG 20 

C 10= INTEGER* 2 ARRAY, 16 LONG, 

C (VALUES 0-4095 , SUBTRACT 2 0 48) *5/20 2 8 . GIVES VOLTAG 

C IR S = NUMBER OF RESINCS ALLOWED (ERRORS) 

C IRE C= COUNTER OF RECORDS (FRAMES OF DATA) 

C BLOCK 512 BITS, 32 BITS = RECORD 

C 800 BPI TAPE UNLABLED 

C IRQ= NUMBER- OF ACTUAL RESINCS (ERRORS) 

C 

c 

INTEGER * 2 10 ( 1 6) , IP ( 1 6 1 
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DATA IRR /O/ 



IF (IREC.EQ.O) IS=0 
I ER= 0 

20 FORMAT (16A2) 

IF (IS.NE.O) 3D TO 50 
READ (IUN,20,END=900) IP 
I REC=IR EC+1 
40 IS=IS+1 

IF (IS. LT.17) 30 TO 50 
READ (IUN,20,END=900) IP 
IS = 1 

IREC = IREC + 1 

50 ICH=IMASK (IP(I3) ,3,0) +1 
C WRITE (6,55) 13 H, IS ,IUN , I REC 

55 FORMAT (' RESYNCING ICH, IS , IUN , IREC ',418) 

C 
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IF (ICH. NE. 1) 33 TO 40 



DO 100 1 = 1,16 
IO{I) = IS HIFT (IP (IS) ,4) 

ICH=I MASK (IP (IS) ,3,0) +1 
IF (ICH. EQ.I) 33 TO 80 
IER=I ER+ 1 

WRITE (6,70) 10 N, IR EC ,I,ICH,IER 

70 FORMAT (' UNIT' , 13 ,’ RECORD 1 , 16 ,' CHAN S DATA CH ',214, 
•ERRORS ',17) 

80 IS=IS+1 

IF (IS. LT.17) 30 TO 100 
READ (IUN,20,EN D=900) IP 
I S= 1 

I REC =IR EC + 1 
100 CONTINUE 
C 
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IF (IER. EQ.O) SO TO 150 



IRR=IRR+ 1 

IF (IRR.LT.IRS) GO TO 120 
WRITE (6,110) 

110 FORMAT ('1 STOPPED IN SUB RD 3ECAUSE 
C OF I RR. GT. ' , 16, ' AT L 1 1 0 * ) 

IRQ= IRR 
STOP 

120 CONTINUE 

WRIT E (6 , 130) I REC, IRR 

130 FORMAT (’RSSYNC AT FRAME ', 16 ,' WITH TOTAL ERRORS’, 17) 
I ER= 0 
IRQ= IRR 
GO TO 50 
150 CONTINUE 
RET U EN 
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900 WRITE (6,910) IUN,IREC 

910 FORMAT ( • 1 END OF UNIT ',13,' AT REC ',17) 
STOP 
END 

FUNCTION ISHIFT (IN,NPLC) 

C RETURNS SHIFTED VALUE OF 1*2 WORD IN 

C -VE LE FT , + VE RIO HT SHIFT 

C 

INTEGER * 2 IN 
IPs IN 

IF (IP.LT.O) IP =IP+ 65536 
IF (NPLC.LT.O) GO TO 30 
ISHIFT=IP/(2**I ABS ( NPLC) I 
RETU RN 

30 ISHIFT=IP*(2**I ABS(NPLC) ) 

IF ( ISHIFT.GT. 6 5535) ISHIFT = MOD (ISHIFT , 65536) 
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RETURN 



END 

FUNCTION IMASK (IN, IBL r IBR) 

C MASK 1*2 WORD IN OUTSIDE BITS IBL S I3R 

C 

INTEGER * 2 IN, 10 
10=1 N 

IF (IBR. EQ. 0) 30 TO 50 
IT=ISHIFT (IN, IBR) 

10= IT 

50 IP=I£HIFT(I0,IBL-15-IBR) 

10= IF 

I MAS K=I SHIFT (10 ,15- IBL) 

RETU 5N 
END 

/* 
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//GO. FT20F0 01 DD UN 11 = 3400 -4 , VDL = S ER=GMDT 1 B, DIS P= (OLD, PASS) , 



// LABEL=(1 ,NL, f IN) , 

// DCB= (RECFM=FB,LRECL=32,BLKSIZE=512,DEN=2) 

//GO. FT21F001 DD UN IT= 3330 V, MS VG P = PUB4 A, DI SP= (OLD , KEEP) , 
// DSN=MSS. S1592.GMDT1A, 

// DCB= (REE FM=VBS,BLKSIZE=4096,LRECL=4092) 

//GO.SYSDUMP DD SYSOUr=A 
//GO.SYSIN DD * 

X COIL PSD, MIRY BA Y, 1 1 22- 1 254 L DC AL , 1 7AUG8 2 , 

20 AVGS, 3 2 S/SEC, 5VOLT , (DB NT**2 VS LOG FREQ) 

Y COIL PSD, MTRY BA Y, 1 1 22- 1 254 L OC AL , 1 7AUG8 2 , 

20 AVGS, 3 2 S/SEC, 5VOLI , (D3 NT**2 VS LOG FREQ) 

STOKES NO. ONE, MTRY BAY, 1122-12 54 
LOCAL, 17AUG82, 20 AVGS, 32S/SEC, 5 VOLT 
STOKES NO. TWO, MTRY BAY, 1122-1254 
LOCAL, 17AUG82, 20 AVGS, 32S/SEC, 5 VOLT 
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STOKES NO. THREE, MTRY BAY, 1122-1254 
LOCAL, 17AUG82, 20 AVGS, 32S/SEC, 5 VOLT 
STOKES NO. ZERO, MIRY BAY, 1122-1254 LOCAL, 
17AOG82, 20 AVGS, 32S/SEC, 5VOLI, (LOG VS LOG) 
PHASE OF X AND Y, MTRY BAY 1 122-1254 
LOCAL, 17AOG82, 20 AVGS, 32S/SEG, 5VOLT 
COHER OF X AND Y COILS, MTRY BAY, 1122-1254 
LOCAL, 17AUG82, 20 AVGS, 32S/SEC, 5VOLT 
X COIL PSD, LA MESA, 11 22-12545L0CAL,17AUG82, 

20 AVGS, 32 S/SEC, 10VOLT, (DB NT**2 VS LOG FREQ) 
Y COIL PSD, LA MESA, 1 1 22-1254 LOCAL, 1 7AOG82 , 

20 AVGS, 32 S/SEC, 10VOL T, (DB NT**2 VS LOG FREQ) 
STOKES NO. ONE, LA MESA, 1 12201254 
LOCAL, 17AUG82, 20 AVGS, 32S/SEC, 10VOLT 
STOKES NO. TWO, LA MESA, 1 122-1254 
LOCAL, 17AUG82, 20 AVGS, 32S/SEC, 10VOLT 
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STOKES NO. THREE, LA MESA, 1 122- 1254 



LOCAL, 17AUG82, 20 AVGS, 32S/SEC, 10VOLT 
STOKES NO. ZERO, LA MESA, 1122-1254 LOCAL, 
17AOG82, 20 AVGS, 32S/SEC, 10VOLT, (LOG VS LOG) 
PHASE OF X AND Y, LA MESA, 1 122- 1254 
LOCAL, 17AUG82, 20 AV3S, 32S/SE3, 10VOLT 
COHER OF X AND Y COILS, LA MESA, 1122-1254 
LOCAL, 17AUG82, 20 AV3S, 32S/SE3, 10VOLT 
RT CIRC POLARIZATION P SD ,MTRY B A Y , 1 1 2 2- 1 254 
LOCAL, 17AUG82, 20 AVGS, 32S/SEC, 5VOLT, (DB VS LOG) 
RT CIRC POLARIZATION PSD, LA MESA, 1122-1254 
LOCAL, 17AUG 82, 20 AVGS, 32S/SEC , 1 3 VOLT , (DB VS LOG) 
LEFT CIRC POLARIZATION PSD, MTRY B A Y , 1 1 22- 1 2 54 
LOCAL, 17AUG 82, 20 A V3S ,32S/SEC , 5 VOLT , ( DB VS LOG) 
LEFT CIRC POLARIZATION PSD, LA MESA, 1122-1254 
LOCAL, 17AUG 82, 20 AVGS, 32S/SEC, 1 0 VOLT, (DB VS LOG) 
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COHER RT CIRC POLARIZATION MTRY BAY/LA MESA, 



1 122-1254 LOCAL, 17AU382, 20 AV3S, 32 S/SEC 
COHER LEFT CIRC POLARIZATION MIRY BAY/LA MESA, 
1122-1254 LOCAL, 17AU3 82, 20 AV3S , 32S/SEC 
/* 

// 
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APPENDIX D 
DATA CURVES 
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Figure D.4 Stokes Parameter One, Mtry Bay, 
17 Aug. 82, 1122-1254 Local 
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Figure D.6 Stokes Parameter Three, Mtry Bay, 
17 Aug. 82, 1122-1254 Local 
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Figure D.7 Phase of x and y Coils, Mtry Bay, 
17 Aug. 82, 1122-1254 Local 
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Figure D.8 Coherence of x and y Coils, Mtry Bay, 
17 Aug. 82, 1122-1254 Local 
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LOG FREQUENCY (HZ) 

Figure D.10 Left Circular Polarization Power Spectral Density 

Mtry Bay, 17 Aug. 82, 1122-1254 Local 
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Figure D.14 Stokes Parameter One, La Mesa 
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Figure D.15 Stokes Parameter Two, La Mesa, 

17 Aug. 82, 1122-1254 Local 
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Figure D.16 Stokes Parameter Three, La Mesa, 
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<r 



Figure D.17 Phase of x and y Coils, La Mesa 

17 Aug. 82, 1122-1254 Local 
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Figure D.18 Coherence of x and y Coils, 

17 Aug. 82, 1122-1254 Local 
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LOG FREQUENCY (HZ) 

Figure D.19 Right Circular Polarization Power Spectral 

Density, La Mesa, 17 Aug. 82, 1122-1254 
Local 
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FREQUENCY (HZ) 

Figure D.21 Coherence of Right Circular Polarization, 

Mtry Bay/La Mesa, 17 Aug. 82, 1122-1254 
Local 
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FREQUENCY (HZ) 

Figure D.22 Coherence of Left Circular Polarization, 

Monterey Bay/La Mesa, 17 Aug. 82, 1122-1254 
Local 
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